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Abstract. We consider hyperfine splitting of Is and, in part, of 2s levels in light hydrogen-like atoms: 
hydrogen, deuterium, tritium, helium-3 ion, muonium and positronium. We discuss present status of pre- 
cision theory and experiment for the hfs intervals. We pay a special attention to a specific difference, 
D21 = 8i?hf s (2s) — _Bhf s (ls), which is known experimentally for hydrogen, deuterium and 3 He + ion. The 
difference is weakly affected by the effects of the nuclear structure and thus may be calculated with a high 
accuracy. We complete a calculation of the fourth order QED contributions to this difference and present 
here new results on corrections due to the nuclear effects. Our theoretical predictions appear to be in a fair- 
agreement with available experimental data. Comparison of the experimental data with our examination 
of D21 allows to test the state-dependent sector of theory of the hfs separation of the Is and 2s levels in 
the light hydrogen-like atoms up to 10 -8 . 

PACS. 12.20.Fv Quantum electrodynamics: Experimental tests - 21.45. +v Few-body systems - 31.30. Jv 
Relativistic and quantum electrodynamic effects in atoms and molecules - 32.10.Fn Fine and hyperfine 
structure 



1 Introduction 

The hyperfine structure (hfs) interval in the ground state 
of a number of simple atoms (hydrogen , deuterium 0] , 
tritium jj| and helium-3 ion Q) has been measured with 
a high precision. The hfs separation in the 2s metastable 
state in hydrogen (5||(|, deuterium and the 3 He + ion 
[^,U was also measured accurately. Some experimental re- 
sults are as old as almost fifty years, but the accuracy of 
even present-day theoretical calculations for the hfs inter- 
val in those light atoms is much lower than that for the 
experiments (see e.g. Table [l]). Theory of the hfs interval 
in simple atoms is essentially based on the bound state 
Quantum Electrodynamics ( QED) , however effects due to 
the nuclear structure are unavoidable and they strongly 
affect the energy levels. Their uncertainty limits the the- 
oretical accuracy for the hyperfine splitting in hydrogen, 
deuterium, tritium and helium-3 ion on a level of 10-200 
ppm. 

One of ways to avoid the problem of the nuclear ef- 
fects is to study atoms free of any nuclear structure such 
as muonium and positronium. Hyperfine split ting in these 
pure leptonic atoms was measured 0,|l^[l2],[F|] with an 
accuracy appropriate for precision tests of the bound state 
QED. Other possibilities to avoid the problem of lack of 



accurate knowledge of the corrections induced by the nu- 
clear effects is related to a fact that those corrections are 
proportional to the squared value of the wave function at 
the origin 



AE(Nucl) = A(Nucl) x \& nl (r = 0)| 5 



■Mr = 0) = 



(Zafm) 



(1) 
(2) 
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where a is the fine structure constant, Z is the nuclear 
charge and m# is the reduced mass of the orbiting parti- 
cle. The relativistic units in which % = c = 1 are used here 
and through the paper. We ignore a difference between en- 
ergy interval E and a measured frequency v = E/h: pre- 
senting the theoretical expression for the energy splitting 
E and numerical results for the frequency v. Here W n i(r) 
is the Schrodinger-Coulomb wave function and A(Nucl) is 
a nuclear parameter which does not depend on the atomic 
state nl. Comparing the hfs for the atoms with a differ- 
ent value of !^ n i(0) one can reduce influence of the nuclear 
structure and test the bound state QED with a high ac- 
curacy. There are two options to vary W n i(0): 

— to compare muonic and electronic atoms (i.e. to study 
atoms with the same nucleus and different values of 
m R ); 
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— to compare hfs intervals for the ns states with different 
value of n or to study the hyperfine splitting for states 
I + 0. 

Presently accurate experimental data are available only 
for one of these two options: it is possible to take advan- 
tage of existence of precision experimental data on the 
Is and 2s hyperfine intervals in a few of light two-body 
atomic systems. A comparison of the data for the Is and 
2s hfs intervals allows to determine value of a specific dif- 
ference 

D 2 i = SE Ms (2s) - E h{s (ls) (3) 

in hydrogen, deuterium and helium-3 ion. The theory of 
this specific difference can be developed much more suc- 
cessfully than that for the ground state interval £ms(1s) 
because of the essential cancellation of the nuclear effects 
(see Eqs. (|, 1). 

The most accurate experimental value for D21 was ob- 
tained for the hclium-3 ion 



D^ P ( 3 He + ) = 1 189.979(71) kHz 



(4) 



after comparison of results obtained for the Is state in 
1969 J! and for the 2s state in 1977 §. The QED theory 
was developed by that time up to third order corrections 
including the (Za) 2 E F , a(Za) 2 E F and (Za) 2 (m/M)E F 
contributions (here Ep is the so-called Fermi energy, lead- 
ing contribution to the Is hfs separation). The experimen- 
tal result in Eq. (||) happened to be in some agreement 
with theory, however, uncertainty of theory was not prop- 
erly estimated. Here we present new theoretical results on 
D21 in hydrogen, deuterium and helium-3 ion ]i4| . In our 
paper we demonstrate that there are a number of higher- 
order QED corrections which were not taken into account 
and which are competitive with the uncertainty of the ex- 
periment. We complete calculation of fourth order correc- 
tions and present theoretical results with accuracy higher 
than that for the measurements. The higher-order nuclear- 
structure effects also contribute to the difference D21 and 
their contribution is important for a comparison with ex- 
periment. They are considered in our paper in detail. 

The paper is organized as following: in Sect. 2 we con- 
sider the QED theory of Is hfs interval and determine 
parameters A(Nucl) for hydrogen, deuterium and 3 He + 
ion. Sect. 3 is devoted to QED calculations of the differ- 
ence Z?2i- We study the fourth order QED contributions 
and, in particular, we find the vacuum polarization con- 
tribution in order a(Za) 3 Ep and the leading logarithmic 
recoil term in order (Za) 3 (m/M)Ep. The nuclear effects 
are taken into account in Sect. 4. We show that study of 
the difference D21 provides an effective test of QED the- 
ory of the hfs intervals -Ehfs(ls) and Ehf s (2s) on a level of 
accuracy essentially below 1 ppm and such a test is free of 
problems of the nuclear structure. That is quite compet- 
itive with investigations of the hyperfine splitting in the 
ground state of muonium and positronium and we present 
a brief overview of them in Sect. 5. Section 6 summarized 
the paper and a comparison of theory and experiment is 
presented there for the difference D21 in hydrogen, deu- 
terium and helium-3 ion and for the ground state hfs sep- 
aration in muonium and positronium. 



2 Hyperfine splitting in the ground state 
in hydrogen, deuterium and helium ion 

The hyperfine splitting of an ns state in a hydrogen-like 
atoms is determined in the non-relativistic approximation 
by the so-called Fermi energy: 



E Ms {ns) = — , 

P /. 8 3 2 a 21+1 ( M 

Ep h = - Z a a c • Ry — — — 

3 /is 21 \m + M 



(5) 
(6) 



Here Ry is Rydberg constant, c is speed of light, h is 
the Planck constant, /is is the Bohr magneton, m is the 
electron mass, M is the nuclear mass and / is the nuclear 
spin. The nuclear magnetic moment (i in our notation can 
be negative (if its direction is opposite to the nuclear spin) 
and the Fermi energy Ep defined as the splitting between 
states with atomic angular moments F = I + 1/2 and / — 
1/2, calculated within the non-relativistic approximation, 
can be negative as well. 

The result of the QED calculations is of the form 



E^ T3 (ls) = E F (l + Q QED (ls)) , 



(7) 



and 



Qwld! US = a. <{ ^{Za) 2 +a{Za) (ln2 - | 



+ 



a(Za) 



— In- -77 fin- 777 

3 (Za) 2 \ (Za) 2 



41n2- 



281 
240 
8, „ 34 
15 ln2+ 225 



17.122 339. 



0.7718(4) 



a 2 (Za) 



.(8) 



The references to all terms can be found in a review in 
Ref. jL8) . The expression above is a result of the external 
field approximation. The recoil corrections involve inte- 
gration over high momentum k ~ M and a consideration 
of a nucleus as a point-like one is not valid in such a case. 
Actually a theory of a point-like particle with an anoma- 
lous magnetic moment is inconsistent and leads to a di- 
vergency at high momentum transfer for the nuclear ver- 
tex. The muon (and electron) must be treated as either a 
point-like particle without anomalous moment, or a parti- 
cle with the anomalous magnetic moment and an internal 
structure. The magnetic anomaly and "structure" effects 
come from the same diagrams. As a result one need to be 
very careful when considering a "point-like" nucleus as an 
approximation. It is important to mention that the lead- 
ing recoil corrections p9[ and most of non-leading terms 
involve the nuclear structure effects and proportional to 
the |*„ ; (0)| 2 (cf. §))• 

To compute any numerical result we use here a 1 — 
137.036 000, c-Ry = 3.289 841 960-10 12 kHz, a e = 1.159 652- 
10~ 3 and parameters of nuclei collected in Table g. The 
values of the fundamental constants and nuclear param- 
eters are based on data taken from (20lEl|, but we keep 
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Table 1. The ground state hfs interval in hydrogen, deuterium, tritium and helium-3 ion. 



Atom 


pexp 


pQED 


TTiCXp riQED 
^hfs _ii hfe 






[kHz] 


[kHz] 


[kHz] 


[ppm] 


Hydrogen 


1420 405.751 768(1), §| 


1420 452 


-46 


-33 


Deuterium 


327 384.352 522(2), § 


327339 


45 


138 


Tritium 


1516 701.470 773(8), [§ 


1 516 760 


-58 


-38 


3 He + ion 


- 8 665 649.867(10), § 


-8 667 569 


1919 


221 



them here only with the accuracy sufficient for our pur- 
poses. The results of calculation are summarized in Ta- 
ble |l|. Leading recoil corrections depend on nuclear struc- 
ture. Some pure QED corrections of higher order are known 
but not included being essentially smaller than uncer- 
tainty of nuclear effects. 

The theoretical calculations above take into account 
only pure QED terms, while the nuclear effects can be 
estimated via a comparison of the experiment and the 
pure QED theory 



A hfs (Nucl) 



£ h E fe P (ls) 



£ Ex P(ls) - gQ ED (lg) 

l<Mo) 2 | 



(9) 
(10) 



The nuclear models or study experimental data on nuclei 
offer another way to find _E Nucl (ls) and ^4(Nucl) and they 
are discussed in part in Sect. ||. 



3 QED calculations of D 2 \ in light atoms 

The evaluation of the QED corrections involves contribu- 
tions of the second, third and fourth order in unit of the 
Fermi energy. The second and third [|^,[l6],[l7j] order 
corrections were calculated some time ago (see Table ||). 
One of the fourth order corrections, (Za) 4 Ep, was also 



found that time 1 26 . Other fourth order terms were found 
only recently and the theoretical expression is now of the 
form 



C 2 i(QED) = {ZafE F x 



5 ^ 111 (7 V 
— (Za) 

8 128 1 ; 



Hi 

y 



ln2-7 ln(Za) - 5.221 23 . 



a 


Am2- 


7 " 






7T 


15 


10 


m 




In 2 


M 


l-H 


~ 




/145 


7 




~ V 128 


~ 8 



1_ 

32 



1 - 



In 2 -7) ln(Za) 



a" ( 16 
2^2 



a 2m ( 16 , „ _ . , 

( T lB2-7]ln(Z«) 



Za m ( 4 , \ , , „ . 

a(Za)(c SE + C VP ) \, 



(11) 



(12) 



where 

/z M 1 
' Hb m Z I 

Two corrections in the fourth order, a(Za) 2 (m/M)Ep 
and a 2 {Za) 2 Ep, were found in Ref. P7J in the leading 
logarithmic approximation and their uncertainties are es- 
timated by a half-value of the leading logarithmic terms. 
The coefficients Cse and Cyp related to the self-energy 
and vacuum polarization higher-order radiative corrections 
were first estimated in Ref. [^7|, but with some misprints. 
Below we correct that estimation and discuss a recent cal- 
culation in Ref. J2£| . 

The expression takes onto account some recoil effects. 
As it was demonstrated by Sternheim ]l6| the n-dependent 
part of the (Za) 2 {m/M)Ep contribution into Ebf s (ns) 
does not depend on the nuclear structure. In the case of 
the a(Za) 2 (m/M)Ep and (Za) 3 (m/M)Ep that is correct 
at least for the logarithmic terms. The pure recoil log- 
arithmic correction (Za) 3 (m/M)\n(Za)Ep is evaluated 
here. A logarithmic part of the QED correction in order 
(Za) 3 m/M is easy to calculate with help of effective po- 
tentials which are responsible for (Za) 5 m 2 /M correction 
to the Lamb shift (cf. f^lll). The result is 



2 (Za) 5 Tfi , 

2 ■ '- In 

3 7T M 



1 

Za 



Ep 



(Zamy 



In 2 



(13) 



However this result is of a reduced value as far as the ef- 
fective potential for the Lamb shift in order (Za) 5 m 2 /M 
has been used. The logarithmic term (with ln(Za)) in that 
order is not dominant. That often happens with pure re- 
coil contributions (in contrast the logarithmic terms dom- 
inates in the case of most of radiative and radiative-recoil 
corrections). The other logarithmic contribution in or- 
der (Za) 5 m 2 /M (with a recoil logarithm \n(mR)) and a 
part of the non-logarithmic term are effectively included 
into the nuclear-structure contributions (see the next sec- 
tion). An essential non- logarithmic part which is not in- 
cluded there is related to two-photon effective potentials 
with derivatives. The scale of the loop integration mo- 
mentum is determined by the electron mass and the re- 
lated contribution does not depend on the nuclear struc- 
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Table 2. Parameters for calculations of the hfs interval in hydrogen, deuterium and helium-3 ion po|pl| . The proton charge 
radius is taken from Ref. 12411 . 



Atom 


Z 


I 


M/m 




E F 


'/ 


Re 










[icr 3 ] 


[kHz] 




[fm] 


Hydrogen 


1 


1/2 


1836.153 


1.521032 2 


1418 840 


5.585 69 


0.88(3) 


Deuterium 


1 


1 


3 670.483 


0.466 975 5 


326 968 


1.714 03 


2.13(1) 


Tritium 


1 


1/2 


5 496.922 


1.622 393 6 


1515 038 


17.831 




3 He+ 


2 


1/2 


5 495.885 


-1.158 750 5 


-8 656 598 


-6.368 36 


1.67(1) 



ture (cf. the Sternheim contribution into Eq. ( |ll| ) in or- 
der (Za) 2 (m/M)Ep). It can be essentially enhanced be- 
cause of a big value of the nuclear anomalous magnetic 
moment and we estimate that non-logarithmic contribu- 
tion as ±r)(Za) 3 /irm/M E F . 



Table 3. QED contributions up to third order to the D21 in 
hydrogen, deuterium and helium-3 ion. 



The higher-order vacuum polarization correction re- 
lated to Cyp is found in this paper. The contribution 
comes from diagrams depicted in Fig. 1. They are eval- 
uated with the exact Dirac wave functions and the Green 
function of an electron in the Coulomb field and the result 
is expanded in powers of Zot. The results for the vacuum 
polarization contribution to the 2s hfs interval are (cf. 
Ref. §§): 



Contribution 


H 


D 


3 He+ 


AE(Fig. la) = 




[kHz] 


[kHz] 


[kHz] 




(ZafE F 


47.222 


10.882 2 


-1 152.439 


+ 


a(Za) 2 E F (SE) 


1.936 


0.446 1 


-37.441 5 




a(Za) 2 E F (VP) 


-0.058 


-0.013 4 


1.4148 




(Zaf%E F 


-0.162 9 


-0.009 4 


-0.796 7 





Total up to 3rd order 48.937 1 11.305 6 -1 189.262 4 



a Ep 
n~8~ 



3tt 



Za 



7 

10 



(Za) 2 



(Za) 



3 A437T 

V 192 



3tt ( Za 



AE(Fig. lb) = 



a Ep 

fT 



3tt 
~8~ 


Za 


34 


8 


225 " 


~ 15 



+ (Za) 



3 / 17157T 7T 



V 1152 



6 ln 



Za 



Finally one can obtain 



+ 



Fig. 1. Vacuum polarization contributions to hfs 



AEH(2s) = ^ x 



(Za) 3 



— Za 

4 




247 
~450 ~ 


A H Za) 


2573tt 


24 V 


1152 



Za 

T 



(14) 



The derivation is considered in detail in Appendix |A|. To 
find a correction to the difference D 2 i one has to compare 
the result for the vacuum polarization contribution to the 
2s hfs obtained above with that for the ground state [E9l 
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5 



Ej? x 



3tt 



Za 



34 



+ < Z °> 225 - 15 lQ(2 ^ 



. „ , q / 5397T 137T , / Za 

+ (ZaY In 

^ 7 V 288 24 V 2 



Finally we find 

At>v p = -E F x { {Za) z ( ~ 



VP a 771 „ J „a2 / 7 8 



10 15 



ln(2) 



(15) 



Ref . |27j . The difference for complete and partial results is 
numerically small for both: the vacuum polarization and 
self energy. That is related to the fact that only the second 
derivative of the wave function at origin depends on n 



(Zam) 



3/2 



7T 1 /2 n 3/2 



x < 1 — (Zamr) + 



(Zamr) z 



1 — n 2 (ZamrY 



6 



(19) 



and the n-dependent coefficient is relatively small. Un- 
der these circumstances we consider the partial results in 
Eqs. ( |l8| ) as a confirmation of direct calculation of the 
self-energy Q and vacuum polarization (see Eq. (p"7|)). 

A summary of the contributions of the fourth order 
terms is presented in Table ||. 



-(Za) 



137T , . . 



(16) Table 4. Fourth order QED contributions to the D21 in hy- 
drogen, deuterium and helium-3 ion. 



and thus 



139 13 



Cvp -ii + 2i ln2 " - 74 - (17) 

As we have mentioned some partial results on a(ZaY E F 
terms were presented in Ref. fl27| with a misprint. The rel- 
ative sign in Eqs. (24,25,27) of paper |27| is to be corrected 
and the result is (the corrected sign is marked by *) 



Csi 



139 „ „ 






41n2 




--ln2 


16 




2 



VP 



5 

24 



13 



-ln2 



-In 2 

* 3 

4 



ln(2) 
In 2 + 



3 

16 



3_ 
16 



Contribution 


H 


D 


3 He+ 




[kHz] 


[kHz] 


[kHz] 


(Za) 4 E F 


0.005 6 


0.0013 


-0.543 


a 2 {Za) 2 E F 


0.003 3(16) 


0.000 8(4) 


-0.069(35) 


a{ZaY%E F 


-0.0031(15) 


-0.000 4(2) 


0.022(11) 


a(ZafE F (SE) 


0.008 3(10) 


0.001 9(2) 


-0.395(37) 


a(ZafE F (VP) 


0.003 


0.000 7 


-0.145 


(ZaY^E F 


0.000 5(5) 


0.000 1 


-0.007(10) 


Total: 4th order 


0.0178(25) 


0.0043(5) 


-1.137(53) 



It was also then expected |^| that the a(ZaYE F results 
found there are likely incomplete. Recently Yerokhin and 
Shabaev directly calculated the self-energy contribution 
psj after our suggestion 

C SE (Z = 1) = 2.07(25) , 
C S e(Z = 2) = 2.01(19) . 

The self-energy result of Ref. [p8| is affected by the higher 
order corrections and thus slightly depends on Z. 

The complete results (17) and (|18] ) indeed disagree 
with the corrected above partial results in Eqs. 



C SE = Z^-I + 51n 2 2^2.5 
64 4 

29 13 
CvP= 64 + 24 ln2 "°- 83 - 



(18) 



and in part it is caused by appearance of effective non- 
relativistic operators with derivatives. Those operators do 
not contribute into logarithmic corrections to ground state 
hyperfine structure [p2j,pl| and were not considered in 



4 Nuclear-structure corrections to D 2 \ 

The leading nuclear-structure corrections to E^ cl (ls) and 
E^ cl (2s), being proportional to the wave function at ori- 
gin (see Eq. @), cancel each other when calculating the 
difference Ei\ - However, some higher-order nuclear effects 
can shift D21 and, in fact, they do. The corrections related 
to the nuclear structure effects can be splitted into three 
terms l27|,pl 



£> 2 i(Nucl) = £>£ + Dg + D c 2 



where 



D A - 



^21 



n c - 

- L/ 21 — 



hr2 + ^)-(ZaY-E^\l S ) 



In 2 • {Zaf{mR E ) 2 E F , 



±-(ZaY(mR E ) 2 E F 



(20) 

(21) 
(22) 

(23) 
(24) 
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Here 



(25) 



is a ratio of quadratic magnetic and electric charge radii. 

Let us discuss origin and accuracy of the nuclear-struc- 
ture corrections. To find the first term (D^) one has to 
somehow determine a value of the nuclear contribution to 
the ground state hfs separation, E^ cl (ls), which contains 
three kinds of terms: 



— nuclear-finite-size effects of order (Za) 3 (mR)EF, where 
R ~ Re ~ Rm', 

— nuclear polarizability corrections; 

— nuclear recoil corrections of order (Za) 3 (m/M) \n(mR). 

The correction for the Is state was studied for hydro- 
gen and deuterium. In the case of the hydrogen the first 
term is dominant and cannot be calculated with accuracy 
better than 20% because of lack of knowledge of the pro- 
ton magnetic form factor at low momentum transfer pi] , 
|23|| . The proton polarizability cannot be successfully esti- 
mated and delivers an essential contribution to the value 
of ££fe cl (ls) in hydrogen. We expect that a theoretical un- 
certainty of nuclear contribution to the Is hfs is at least 
20% of its value. 

Deutron is a weekly bound nucleus and the dominant 
nuclear effect for the hyperfine separation in the ground 
state of deuterium is related to the deutron polarizabil- 
ity. The nuclear correction was estimated in Ref. |Q as 
43 kHz, however, the uncertainty is not presented there. 
We expect that the uncertainty lies between 10 and 30 
kHz. Our assumption is based on examination of the log- 
arithmic approximation used in Ref. Let us concen- 
trate our analysis on two corrections of —19 and +11 kHz 
which were found in the logarithmic approximation. They 
are proportional to \n(m p / k), where n ~ 45.7 MeV is the 
inverse deutron size and m p is the proton mass. The va- 
lidity of the logarithmic approximation suggests that the 
logarithm is big enough, but that is not really a case: 
ln(m p /«;) ~ 3.0. We expect that the uncertainty of such 
an approximation lies between 10 and 30 kHz, that de- 
pends on possible correlations between these two logarith- 
mic contributions. 

To the best of our knowledge there are no results pub- 
lished on the nuclear contributions to the hyperfine sepa- 
ration in the tritium atom and the helium-3 ion. Due to 
lack of accurate calculations for i^ f " cl (ls) we estimate the 
nuclear structure contribution to the Is hfs interval in all 
atoms discussed above by a comparison of experimental 
data with a result of the QED calculations (see Table [j] 
and Eqs. (§ | i ©■ 

All these corrections to -Ehf" cl (ls) are related in the 
leading order to the two-photon exchange with a hard- 
momentum exchange loop (fc 3> Zam). Their calculation 
is similar to that for a(Za) 3 Ep terms. Some of these two- 
photon contributions can induce additional terms with 
derivatives but that will involve an additional suppressing 
factor m/k. The factor m/k is small for the finite-size and 
polarizability and for a logarithmic part (with ln(mi?)) 
of the nuclear recoil contribution. It is about unity only 



for a part of the nuclear recoil contribution related to low 
momentum transfer k ~ m and k < m, however, we have 
calculated those corrections of order (Za) 3 (m/M) in the 
leading logarithmic approximation in Eq. ( |l3|) . We expect 
these contributions are relatively small, because the recoil 
effects are not dominant in the term and because of 
small numerical coefficient for state-dependent terms in 
the hydrogenic wave function at origin (see Eq. ( |l9| ) ) . The 
7-matrix structure is close to that in the case of the vac- 
uum polarization where the terms with derivatives induce 
numerically small contributions. We estimate the uncer- 
tainty of such an approximation for _D^\ as 10%. 

To verify the expression for we also compare value 
' Is) with that for E^ cl (2s)/8 (see Table |). The 
latter was found via a comparison of a pure QED theoret- 
ical expression (cf. Eqs. (|^,|Tl|) 

E^ D (2s) = ^f(l + Q QED (2s)) . 



oiE^ 



Qqed(2s) 



^(Za) 2 + a(Za)[\u2- - 



a(Za) 



81n2- 

247 
~450 



_2 1 

3 11 {Zaf 



In - 



1 



1541 
240 



(Za) 2 
11.901106... 



0.7718(4) ■ 



\Za)\ 



(26) 



withexperimental data. The results for the nuclear contri- 
butions to the Is state in Table [l] and for the 2sstate in 
Table | agree with each other. 



Table 5. 2s hyperfine splitting in light atoms. Zi-Ehfs = 
£§f D ) and measured in ppm in respect to Ef/S. 



Atom, 


piCXp 

^hfs 


tttQED 
^hfs 


AE Ma 


state 


[kHz] 


[kHz] 


[ppm] 


H, 2s 


177 556.785(29), § 


177 562.7 


-33 


H, 2s 


177 556.860(50), || 




-32 


D, 2s 


40 924.439(20), § 


40 918.81 


137 


3 He+, 2s 


-1083 354.981(9), § 


-1083 594.7 


221 


3 He+, 2s 


-1083 354.99(20), § 




221 



Two other nuclear contributions, and DJ?i (see 
Eq. ( p2| ) and Eq. are smaller than and their 

evaluation is similar to that for the a(Za) 2 EF contribu- 
tions and completely understood. They were derived in 
Refs.^7],^3| with help of some effective potentials, and the 
result does not depend on any nuclear models. The result 
for the nuclear contribution in Eq. (|20l) can be presented 
in a form slightly different from Eqs. (|20|, |2l| , ^g , |23|) 



(Zar-AE^\ls) 



£> 21 (Nucl)= (ln2+ — 
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— - 21n2 
8 



^Lamb 



(Is) 



-5 C 

8 s (Za) 2 m 



(Za) 2 m 

^F j 



Table 8. Muonium hyperfine structure 



which is more useful for phenomenological applications. 
We check this expression for the effects caused by a distri- 
bution of the nuclear charge and magnetic moment within 
some models in Appendix [B| and confirm it. The results on 
nuclear contributions to the D21 in the light hydrogen-like 
atoms are presented in Table a 



Table 6. Nuclear-structure contributions to the D21 in hydro- 
gen, deuterium and helium-3 ion. 



Contribution 


H 


D 


3 He+ 




[kHz] 


[kHz] 


[kHz] 


D& 


-0.002 2(2) 


0.002 1(2) 


0.360(36) 




0.000 3 


0.000 4 


-0.028 5 


Dg 


-1 ■ 10 -4 • c 


-1.3-10 _4 -C 


8.6 ■ 10 -3 ■ C 


D 2 i(Nucl) 


-0.002 


0.002 6(2) 


0.332(36) 




-1 ■ 10~ 4 ■ c 


-1 ■ 10~ 4 ■ c 


+9 • 10" 3 ■ C 



The final theoretical results 

D 21 (theor) = £> 21 (QED) + D 21 (Nucl) 



(27) 



for hydrogen, deuterium and helium-3 ion are summarized 
in Table 0. The table contains also the experimental re- 
sults. The main sources of uncertainty of the theoretical 
calculations are related to 

— the use of logarithmic approximation in evaluation of 
higher order QED corrections of order a 2 (Za) 3 E F , 
a(Za) 2 (m/M)E F and (Za) 3 (m/M)E F ; 

— calculation of higher-order nuclear effects. 

The recoil contributions in order a(Za) 2 (m/M)E F and 
(Za) 3 (m/ M)E F also limit accuracy of the calculations of 
the hyperfine splitting in the ground state of muonium and 
positronium and we overview theory of these two quanti- 
ties in the next section, 
table 6 



Term 


Fractional 


AE 




contribution 


[kHz] 


Fermi energy 


1.000 000 000 


4.459 031.92(51) 




0.001 159 652 


5 170.93 


2nd order QED 


- 0.000 195 815 


- 873.15 


3rd order QED 


- 0.000 005 923 


- 26.41 


4th order QED 


- 0.000 000123(49) 


- 0.55(22) 


Hadronic effects 


0.000 000 054(1) 


0.24 


Weak int. 


- 0.000 000 015 


0.06 


Total 


1.000 957 830(49) 


4463 302.91(51) (22) 



7,32 



J ) and we follow consideration 



J34| (see also Refs. 

there. The hadronic contribution is taken from Ref. 
|35| . The results are collected in Table ||. 

The dominant QED contribution to the uncertainty 
(0.22 kHz) comes from unknown non-leading terms in or- 
ders a(Za) 2 m/M and (Za) 3 m/M, which are estimated 
by a half-value of the leading double logarithmic correc- 
tions |32,31|, despite some terms beyond the double loga- 
rithms are known (see Ref. |34| for discussion). 

To find an absolute value one has to determine the 
Fermi energy in Eq. @ which contains the magnetic mo- 
ment of the muon. The most accurate value of it can be 
obtained from study of the ground state hyperfine struc- 
ture in the magnetic field |M. The related uncertainty is 
0.51 kHz. Opposite to the theory of -D21, the Fermi energy 
has to be calculated very accurately and its value depends 
on our choice of a value of the fine structure constant. Here 
we use an original result from study of anomalous mag- 
netic moment of the electron a" 1 = 137.035 999 58(52) 
[p8| . The related uncertainty is only 0.03 kHz, however, 
scattering of various results for the fine structure con- 
stants (see e.g. Ref. |2C|0 corresponds to a much bigger 
uncertainty. 

The positronium hfs interval can be calculated and 
measured less accurately than that in muonium. However, 
it provides us with a sensitive test of the same recoil cor- 
rections as in the case of muonium and -D 2 i- The recoil 
effects in positronium are essentially bigger than in other 
atoms, because M = m, and they may be studied in de- 
tail. The positronium hyperfine splitting in the ground 



5 Hyperfine structure in pure leptonic atoms state can be presented in the form 



The theoretical expression for the hfs interval in the muo- 
nium ground state can be presented in the form 



£ h fs(theor) = E F (1 
= E F {1 



Q) 



Q2 + Q 3 + Q4 + Qh + Qi, 



and the results for the QED contributions of the second 
(Q2), third (Q3) and fourth order (Q4), for the hadronic 
(Qh) and weak contributions (Q w ) are reviewed in Ref. 



.Ehfs(theor) = E F (l + q 1 a + q 2 a 2 + q 3 a 3 ) 
= E F (1 + Qr + Q 2 + Q 3 ) , 

where coefficients slightly depend on a containing In a. 
The Fermi energy is defined in positronium in a different 
way (comparing with Eq. (^|) for hydrogen and muonium) 



7 

E F (Ps)/h = -a 2 cRy 
6 



(28) 
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Table 7. Value of D21 in hydrogen, deuterium and helium-3 ion. Results for nuclear correction and theory are for f = 0. 
Z\(exp — th) = D2i(exp) — D2i(theor) and a is a final uncertainty of Zi(exp — th). The final nuclear contribution D2i(Nucl) 
here is presented at f = 



Value 


H 


D 


3 He+ 


D 2 i(exp) [kHz] 


48.53(23), @ 


11.16(16), § 


-1 189.979(71), § 


Au(QED) [kHz] 


48.955(3) 


11.309 9(5) 


-1 190.400(53) 


Au(Nucl) [kHz] 


-0.002 


0.002 6(2) 


0.332(36) 


D 2 i(theor) [kHz] 


48.953(3) 


11.312 5(5) 


-1 190.067(64) 


Z\(exp - th) [kHz] 


-0.42(23) 


-0.15(16) 


0.09(10) 


Z\(exp - th)/cr 


-1.8 


-1.0 


0.9 


a/E F [ppm] 


0.f6 


0.49 


0.01 



because of annihilation effects and a symmetric treatment 
of magnetic moments of the electron and the nucleus (posi- 
tron) . The results are summarized in Table ^| (see p9|,[40|, 
[|l], [42, 18, 3(J and references there). The third order correc- 
tions appear to be large because of a double logarithmic 



enhancement [pll (In 



24). The value of q% is calcu- 



lated with taking into account recent results on a' 'm In a 
correction |36| , however, the uncertainty is estimated by 
a half- value of the leading a 7 m In 2 a H (see Ref. @ for 
discussion). 



Table 9. Positronium hyperfine splitting 



Term 


q 


Q 


AE [MHz] 


Fermi energy 


1 


1.000 0000 


204 386.6 


1st order QED 


-0.674 16 


- 0.004 919 6 


-1005.5 


2nd order QED 


1.084 


0.000 057 7 


11.8 


3rd order QED 


-15.6 


-0.000 0061(22) 


-1.2(5) 


Total 




0.995 132 1(22) 


203 391.7(5) 



6 Summary 

The results for the precision calculations of the hyperfine 
structure in the light hydrogen-like atoms are summarized 
in Table One can see that investigations of the differ- 
ence D21 provide very accurate tests of the bound state 
QED calculations. We consider study of D 2 i as a test of a 
state-dependent sector of theory of the hyperfine splitting 
of the Is and 2s states and so the fractional accuracy of 
such theory is related to the Ep, the leading contribution 
to the Is hfs. The accuracy of comparison of theory and 
experiment can be characterized by a standard deviation 
a which contains contributions to uncertainty from both: 
theory and experiment. The final uncertainty is found to 
be for D21 as small as few part of 10 7 in the case of hydro- 
gen and deuterium and even better in the case of helium 



ion: a part of 10 s . That is competitive with other tests 
of the bound state QED and in order to clarify advan- 
tages and disadvantages of studying D21 let us list main 
problems which theoretical calculations have met by now: 

— there are two essential problems of the bound state 
QED: 

— evaluation of higher-order recoil corrections (that 
is mainly a problem of all QED calculations for the 
hyperfine structure including -D21); 

— evaluation of higher-order two-loop corrections (that 
is rather a problem of the Lamb shift calculation 
and only one value related to the hyperfine struc- 
ture, D21, is sensitive to such corrections); 

— there are two other problems related to other part of 
physics: 

— determination of the fundamental constants (like 
e.g. determination of the fine structure constant 
and magnetic moment of muon needed to calculate 
the Fermi energy Ep); 

— nuclear structure, which affect energy levels and, in 
particular, shifts values of the Lamb shift and the 
hyperfine separation. 

The difference -D21 happens to be an only value that 
is sensitive to both higher-order corrections: recoil and 
two-loop and that is not sensitive to problems beyond 
QED (determination of the fundamental constants and 
nuclear structure). Tests of QED are sometimes consid- 
ered as a search for new physics beyond the Standard 
model. Such exotic contributions are rather expected to 
be proportional to 1/n 3 and must vanish for D21 in the 
leading order. The next-to- leading terms could contribute 
but they effectively would be taken into account in D 2 sr 1 ucl 
being included into i?^ f " cl (ls). That fact makes the differ- 
ence -D21 useful for a very specific test of the bound state 
QED, a test which involves no problem beyond QED. 

Our theoretical predictions appear to be in a fair agree- 
ment with four of five accurate measurements (see Ta- 
bic |Io| ), while a minor discrepancy of 1.8 a with the most 
recent result from Ref. Q is observed. Because of agree- 
ment with other data and espetially with the most accu- 
rate result for helium ion [|| we expect that the problem 
of this minor discrepancy comes from the experimental 
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Table 10. Hyperfine splitting: precision tests of the bound state QED. The final uncertainty a includes contributions from 
both: theory and experiment. References for the D%\ are presented for the both states: 2s and Is. 



Atom 


Value 


Exp. 


Theor 


A 


A/a 


a/E F 






[kHz] 


[kHz] 


[kHz] 




[ppm] 


Mu 


S hfe (ls) 


4 463 302.78(5), 0] 


4463 302.91(56) 






0.12 


Ps 




203 389.1(7) • 10 3 , Q 


203 391.9(5) ■ 10 3 


-2.8(9) 


-3.3 


3.4 


Ps 


£hfe(ls) 


203 387.5(16) ■ 10 3 , Jl|] 




-4.4(17) 


-2.6 


7.9 


H 


-D21 


48.53(23), §/@ 


48.953(3) 


-0.42(23) 


-1.8 


0.16 


H 


D21 


49.13(40), |/@ 




0.18(40) 


0.4 


0.28 


D 


D21 


11.16(16), §/g 


11.312 5(5) 


-0.15(16) 


-1.0 


0.49 


3 He+ 


D21 


-1 189.979(71), §/§ 


-1 190.068(64) 


0.09(10) 


0.9 


0.01 


3 He+ 


£>21 


-1190.1(16), §/§ 




0.03(160) 


-0.02 


0.18 



side. One can see that there have been no improvement 
in microwave measurements of the 2s hfs for the last few 
decades. We expect that some progress is still possible 
and that it is now also possible to perform an optical 
measurement of this quantity via comparison of different 
Is — 2s transitions in hydrogen and deuterium, some of 
which were measured recently very precisely J43J . 
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A The vacuum polarization contribution to 
the 2s hyperfine splitting 

An exact relativistic expression for the vacuum polariza- 
tion correction to hfs in the ground state of a hydro gen -like 
atom with a point-like nucleus was derived in Ref. |29|,|30|. 
Using the same method we can calculate the vacuum po- 
larization contribution for the 2s-state. As in the case of 
the Is state [ p9| , we study a more general case consider- 
ing an orbiting particle with the mass m different from 
the electron mass m e which is related here to a particle in 
the vacuum loop (see. Fig. 1). That offers an opportunity 
to perform some additional tests of our results. 

The diagrams contributing to the hfs separation are 
presented in Fig. 1. We obtain 

a 



2 - 5e + 2e 2 - 2£ 2s {\ + £ 2s ){h - 2e) 
x |-4£ 2s (l + £ 2s )(3 - 2e) 2 (l - e) Jio(^) 



- (1 + 6£ 2s + 6£f J(3 - 2e) 2 (l - 2e) J 20 (k 2 ) 
+ 2(1 + 2£ 2s ) 2 (l - e)(5 - 13e + 6e 2 ) J 30 (n 2 ) 

- (1 + 2£ 2s ) 2 (l - e)(3 - 8e + 4e 2 ) J 40 (k 2 ) 



for the single-potential contribution related to Fig. la and 

o 



Z\£ 2 (VP-U-T) 



E T {2s) 



£ 2s (Za) 2 (3-2e) 2 



e(l-2e) 2 (5 + 4£ 2s -2e)(2-e) 2 
2(l + e)(2 + 2£ 2s -5e + 2e 2 ) 



4 - 9e 2 + 4e 3 
1 - e 



1 - 2e 



J2o(k 2 ) 



^10(^2) 



4-21e + 23e 2 



(1-e) 2 
2(-18 + 31e-19e 2 + 4e 3 ) 
3-2e 
2(-6 + 9e- 10e 2 + 4e 3 



£is ^20(^2) 



^30(^2) 

<?2s <^3o(^2) 



- 2 



3 - 2e 

10 - 9e + 2e 2 + £ 2s (13 - 14e + 4e 2 ) J 40 (k 2 ) 



+ 2(l + 2£ 2s )(2-ef J 50 (k 2 ) 



4(2-e)(l + £ 2s -e) 



^21(^2) 



AE 2 (VP-TU) = -E T {2s) 



1 - e 

- 8(l + £ 2s )(2-e) J S1 (R 2 ) 
+ 4(l + 2£ 2s )(2-e) J A1 {k 2 ) 



for the double-potential term in Fig. lb. Here we mainly 
follow notations of Refs. H,P9| and, in particular, we in- 
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traduce the relativistic Fermi-Breit energy 



E T (2s) = E F x 



2{Zaf 



(l+£ 2s )(5-2e)-l 



(1 + 2£ 2a ){l - e)(2 - e)(3 - 8e + 4e 2 ) 



where 



l_Vl-(Za) a ) 



(Za) 2 



(Za) 



2 

Za to 



k 2 = 



2£ 2 .« 



-2 M 



(Za) 



(Za) 2 



Basic integrals J mn arc defined as 
v 2 (l-v 2 /3) I 



J mn {K>) = dv 



1-v 2 



K\/\ — V 2 

1 + nyl — v 



x In" 



1 + kvT 



They can be expressed in terms of the beta-function and 
the hypergeometric function as follows |44|,p9t 

4o = \n m 5(3/2, m/2) 
x 3 F 2 (m/2, to/2 + 1/2, m/2; 1/2, to/2 + 3/2; k 2 ) 

- y K m+1 B(3/2,m/2 + l/2) 

x 3 F 2 (to/2 + 1, to/2 + 1/2, to/2 + 1/2; 3/2, to/2 + 2; k 2 ) 

- i K m B(5/2,TO/2) 

x 3 F 2 (to/2, to/2 + 1/2, to/2; 1/2, m/2 + 5/2; k 2 ) 
+ 2! k ™+ 1 j B(5/2,to/2 + 1/2) 

x 3 F 2 (m/2 + l, to/2+1/2, to/2 + 1/2; 3/2, to/2 + 3; k 2 ) 
and 



^ m.rj. 



d n J m 



In the case of an electronic atom (to = m e ) and small 
Za we arrive to a result in Eq. (|T^), which reproduces 
all known terms of the expansion over Za E^,^5[|3^] and 
presents a new contribution in order a(ZayEF- 

One can also study the vacuum-polarization contribu- 
tion to the hfs in muonic atoms putting to = m„. In the 
case of low Z and arbitrary n = Zam,, /m r , ~ 1.5Z we 
reproduce the non-relativistic limits J|6). For the case of 



small Za and large k the result can be presented as an 
expansion over Za and k . 



AE 2 (VP) = ^ 

IT O 



+ (Za) 2 



8, , . 4tt 2 85 49 
3 ln(K) + — "18 + ^ 
17, . . 3tt 2 37 

y ln(K) + — -y 



3 4«r 



where ip( z ) is the logarithmic derivative of the _T- function. 
The logarithmic part of the correction can be easily found 
within the effective charge approach (cf. Q,^9|), results 
agree with ( p9|) and that is an additional confirmation of 
our result. 



B A model-dependent calculation of the 
finite-nuclear-size corrections to energy levels 



Here we study the contributions to E^ h (ls), E^ c \ls) 
and D21 related to the distribution of the nuclear charge 
and the nuclear magnetic moment. The distribution is de- 
scribed by the nuclear electric and magnetic form factors 



G a E (q 2 ) = 1 + t 



A 2 



A 2 



1 



Parameters A and p are free parameters, however, we con- 
sider here only linear in (3 contributions. Since we intend 
to verify a model independent expression in Eq. (|27]) , it is 
not important that the distribution above is not quite a 
real one. 

We perform calculations with the Dirac wave functions 
and expand the results over (Za) and m/A. In the case of 
a = 1 the nuclear corrections are 



TjiNucl 
^Lamb 



(Is) 



4(Za) 4 /3 B 



m 



<s Cl (l*) 



£>2i(Nucl) 



-4Za 
3 

4 ' 



m 

(3 E m f} M m 



A E 
41n2 



A 



M 



(Za) 



E f 



Pettl Petti 



A, 



-(12-81n2)(Za) 2 /? £ 



3 , s0 /to 
2 (Za) ^ (aZ 



A m 



E F 



E F 



and 



Re 



Rm 



Q(3e 
A 2 

ly E 

6/3. 



M 



A 2 
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We note that in the case of Pe — Pm — lj Ae = Am we 
arrive at a — 2 to the so-called dipole model commonly 
used as an approximation for the proton internal struc- 
ture. However, since we calculate the linear in j3 terms the 
well known result for the hydrogen hfs cannot be repro- 
duced. In the case of a — 2 we find 



Petti f3 M m 
A E A M 



E^\ls) = -6Za 



9 



£> 21 (Nucl) = - I £ + 6 In 2 ) {Zaf [ + | E f 



Ae Am 

2 

m 



-(24- l&]xi2){ZaYP E ( jj- ) Ep 



-3(Za) 2 p M [ — ] E, 



and 



Re 



Rm 



12fe 

A\ 

YIP. 



M 



A 2 



The results for a = 1 and a = 2 confirm expression 
in Eq. ( p7| ) . Similar calculations can be performed for any 
integer value of a. Let us also mention, that in principle 
any moments of a real distribution ((R n pj^(R))) can be 
reproduced by a finite sum of G E (q 2 ) and G M (q 2 ) over 
integer a after adjusting parameters A(a) and P{a). 
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